
A METHOD OF DETERMINING MODAL DATA 
OF A NONUNIFORM BEAM WITH EFFECTS 
OF SHEAR DEFORMATION AND ROTARY INERTIA 

by Vernon L. Alley, Jr., Robert J. Guillotte, 

and Lessie D. Hunter •si’f/’N 

Langley Research Center 
Langley Station, Hampton, Va. 

NATIONAL AERONAUTICS AND SPACE ADMINISTRATION . WASHINGTON, D. C. • SEPTEMBER 1965 



I 


TECH library kafb, nm 

I linn III.. .... s 



TECH LIBRARY KAFB, NM 



□1547b? 


NASA TN D-2930 


A METHOD OF DETERMINING MODAL DATA OF A NONUNIFORM BEAM 

WITH EFFECTS OF SHEAR DEFORMATION AND ROTARY INERTIA 

By Vernon L. Alley, Jr., Robert J. Guillotte, 
and Lessie D. Hunter 

Langley Research Center 
Langley Station, Hampton, Va. 


NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 


For sale by the Clearinghouse for Federal Scientific and Technical Information 
Springfield, Virginia 22151 — Price $2.00 



A METHOD OF DETERMINING MODAL DATA OF A NONUNIFORM BEAM 


WITH EFFECTS OF SHEAR DEFORMATION AND ROTARY INERTIA 

By Vernon L. Alley, Jr., Robert J. Guillotte, 
and Lessie D. Hunter 
Langley Research Center 


SUMMARY 


A recurrence solution including the secondary effects of rotary inertia 
and shear deformation on discontinuous nonuniform beamlike structures is pre- 
sented for obtaining highly descriptive free-free natural mode data. Results 
of other studies are included to ascertain the significance of the secondary 
influences on the classical uniform beam. Numerical results are also included 
for the application of the method to a first-stage and fourth-stage configura- 
tion of a typical space research launch vehicle. The numerical results indicate 
that shear deformation is generally the prime contributor of the secondary 
effects on typical launch vehicle configurations ; the major influences of the 
secondary effects are seen in reductions in natural frequencies and in changes 
to the mode slopes, mode moments, and' mode shears. 


INTRODUCTION 


Natural vibration characteristics are frequently required for systems that 
can be appropriately represented as free-free beamlike structures. Such struc- 
tures are consistently encountered in the array of stages that comprise the 
typical launch vehicle for space research. 

An assessment of the natural vibration characteristics normally will 
include the calculations of the fundamental frequencies of oscillation along 
with a number of overtones. The deflected curves associated with each char- 
acteristic frequency, known as either mode shapes, characteristic functions, or 
eigenfunctions, are also desired data. In addition, the slopes, moments, and 
sometimes the shears associated with the mode shape are provided. Special 
definite integrals also are frequently computed, for example, the generalized 
or effective mass of a mode. 

These data have a variety of uses. Their most familiar use is in providing 
knowledge of natural frequencies for designing in order to avoid states of 
resonant vibrations. This application is particularly important in designing 
or qualifying the spin program for unguided launch vehicles. (Trajectory dis- 
persion control is frequently accomplished by spinning a vehicle about its 
longitudinal axis.) High dynamic stresses and even structural failures might 



result if the roll frequency is near or coincident for an appreciable time to 
one of the natural frequencies of the system. 

In addition, the vibration response levels of a vehicle will generally be 
high near the frequencies of the natural vibrations of the basic structure. 
Knowledge of the probable distribution of the disturbance spectrum can provide 
information for appropriate shock mounting of delicate payloads and can lead to 
a proper choice of the dynamic characteristics of instrumentation. 

Furthermore, calculated mode shapes and frequencies can be useful in 
interpreting measured vibration data. For example, moment and shear data asso- 
ciated with calculated modal responses are frequently used for rapid assessment 
of the load significance of recorded flight data. 

Probably the greatest value of modal data is realized in series solutions 
of the differential equations of motion of the system for which they have been 
generated. The well-known orthogonality properties of mode shapes, the fact 
that the modes satisfy the appropriate boundary conditions for their particular 
constraints, and the fact that the usual response of a beam system is adequately 
described by superposition of a few modes make them ideal functions for such 
applications, widely known as modal form solutions. 

The modal form series solution is particularly adaptable to formulating 
techniques for yielding gust, wind loads, ignition and separation responses, 
and other transient loads . The modal form approach has also proven popular for 
developing the characteristic equation of structures coupled with closed-loop 
autopilot systems . Such analyses permit investigations of the stability bound- 
aries of such systems and lead to proper gain levels and filter characteristics 
to yield stable performances. 

Such valuable applications of modal data have stimulated analysts to 
develop a variety of analytical techniques to achieve increased accuracy and 
scope of output with a minimization of input effort. Well-known and fundamental 
considerations of beam vibration are presented by Den Hartog in reference 1. 

The classical techniques of Rayleigh and Stodola are two of the methods dis- 
cussed in detail by Hartog. In reference 2, Scanlan and Rosenbaum describe the 
classical procedure attributed to Rayleigh and Ritz along with the Myklestad 
method. They proceed to meticulously outline a step-by-step matrix method 
employing influence coefficients. This latter procedure has been formulated 
into a rigorous matrix solution and is presented by Alley and Gerringer in ref- 
erence 3 . Houbolt and Anderson have organized the method of Stodola and pre- 
sented other useful related material in reference 4. An interesting integral 
series solution technique that is also an effective solution to the nonuniform 
beam vibration problem is presented by Spector in reference 5« These references 
are but a few of the publicized methods relating to natural vibrations of non- 
uniform beams. 

When the analyst is confronted with a particular assignment for computing 
modal data, the choice of the most suitable method should be made in view of the 
objectives of the assignment and the virtues and limitations that characterize 
the methods. Energy methods following the concepts of Rayleigh or the 
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Rayleigh -Ritz approach generally yield adequate results on frequencies hut 
possess weaknesses in the exactness and descriptiveness of the mode shapes and 
their derivatives. The classical Myklestad method as well as the matrix method 
of reference 3 treat the system as a number of lumped point masses coupled with 
massless springs. These discrete mass models also yield adequate frequency 
data but suffer a loss in descriptiveness of the other modal data, particularly 
in the shear curves associated with the modes. The Stodola process, or one 
similar, requires iteration to converge on both frequency and mode shape. 

Whereas the frequency is readily obtained to a desired accuracy, the mode shape 
is not so readily made everywhere convergent to a prescribed accuracy. In addi- 
tion, the procedure requires sweeping techniques to obtain modes of higher fre- 
quencies than the fundamental. This operation frequently presents practical 
problems due to degeneration in numerical accuracy. The integral series method, 
although capable of yielding both the fundamental and higher modes, presents a 
formidable computing problem, the maintenance of sufficient accuracy in per- 
forming numerically the required multi -integrations . In the original or classi- 
cal form, all of the aforementioned methods omit the secondary influences of 
rotary inertia and shear deformation. Many extensions to the classical forms 
have been made throughout industry and government, yet few satisfactory methods 
incorporating the secondary influences have been adequately described in gen- 
erally available literature. Shear deformation introduces an additional source 
of deflection to the customary flexural deformation considered in elementary 
beam theory. Cross-section rotary inertia provides additional dynamic loading 
to the system due to the rotational acceleration of the cross section of the 
beam. In most slender beams the shear and rotary inertia contributions to 
loading and deflection are small in comparison with those resulting directly 
from bending. In multistage launch vehicles of conventional fabrication the 
shear deformation is secondary to bending and the structure is generally so 
slender that rotary inertia may also be ignored. 

However, the use of fiber-glass stages has given unexpected emphasis to 
the importance of shear. The low ratio of shear modulus to the flexural modulus 
of elasticity of fiber glass increases the relative contribution of shear defor- 
mation to bending and thus produces significant shear effects in geometries 
that otherwise could be investigated adequately without considering secondary 
effects. 

In addition, when computing the modal data of beam systems with low length - 
to-diameter ratios such as characterized by the upper stages of conventional 
multistage launch vehicles, the inclusion of both shear deformation and rotary 
inertia has been found to be important to data accuracies. 

This paper presents another technique for computing modal data on nonuni- 
form beams. The method has proven very satisfactory for accurately describing 
not only the mode shape, but its slope, mode moments, and mode shears. The 
method is inherently applicable to structures exhibiting numerous discontinu- 
ities in their mass and stiffness properties and requires no discreting to an 
analogous lumped mass system. The differential equations of motion of the beam 
system are dealt with directly. Modes need not be computed consecutively as is 
required by a number of other methods. The technique has been employed suc- 
cessfully in obtaining accurate modal data up to the tenth mode on complex 
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structures. Furthermore, the formulation accurately incorporates the secondary 
influences of rotary inertia and shear deformation. 

The organization of the problem in the matrix form that follows has proven 
to he highly practical in obtaining numerical solutions on a high-speed digital 
computer. Also, the numerical integration technique that is an integral part 
of the method has been found to be highly stable and accurate. 

The detailed mathematical developments of the method are presented with 
information pertinent to practical numerical solutions utilizing the digital 
computer. Some data are supplied for assessing the probable importance of the 
secondary influences. An example of the application of the method to a typical 
four-stage launch vehicle is submitted with a study of the importance of sec- 
ondary influences on the upper stage. 


SYMBOLS 

A cross-sectional area in shear, in.^ 

"A" solution for superposition method 


Al 

value of 


at 

x z 

due to a unit value of 

* 

at 

x a, in. -1 

a 2 

value of 

V 

at 

X Z 

due to a unit value of 

♦ 

at 

x a , lb/ rad 

a 

coefficient 

(eqs 

. (30)), lb-rad^/in 

2. 



"B" 

solution 

for 

• superposition method 




®1 

value of 

t' 

at 

X Z 

due to a unit value of 

5 

at 

x a , rad/in.^ 

b 2 

value of 

V 

at 

X Z 

due to a unit value of 

5 

at 

x a , lb/in. 

b 

coefficient 

(eqs, 

• (30)), Z n+ ]_a£, lb-rad^/in, 

,2 



C 1 

value of 

£ 

at 

X Z 

due to a unit value of 

+ 

at 

x a , in . / rad 

C 2 

value of 


at 

X Z 

due to a unit value of 


at 

x a , unities s 

D 

maximum diameter 

of beam, in.; also used as determinant (eq. ( 32 )) 

D 1 

value of 


at 

x z 

due to a unit value of 

c 

at 

x a , unities s 

d 2 

value of 

* 

at 

x z 

due to a unit value of 

5 

at 

x a , rad/in. 
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d 

E 

e 

f(co) 

G 

I 

K 

k r 

k 

L 

M 

m 

m r 


p n+l 

P 

t 

V 

X 

x 

x a 

X Z 

Ax 

Z 


shear deformation coefficient (eqs. ( 30 )), lh -1 

modulus of elasticity in bending, lb/in. 2 

flexibility coefficient (eqs. ( 30 ))j -gp Xb~^-ln.~2 

frequency equation (eq. (37)) 
modulus of shear, lb /in. 2 

area moment of inertia for a given cross section, in.^* 
coefficient for cross section in shear 

effective spring constant of the rth mode; Ib/in. when mode 
shapes are considered dimensionless 

half increment (eqs. ( 30 )), Ax/2, in. 

overall length of beam, in. 

mode moment due to bending, in-lb 

mass per inch of length, lb-sec 2 /in. 2 

effective mass of the rth mode; lb-sec 2 /in. when mode shapes 
are considered dimensionless 

matrix in recurrence solution of equation ( 27 ) 

point in derivation 
time, sec 

mode shear due to bending, lb/ in. 
column matrix of variables (see eq. (l0)) 
coordinate along length of vehicle, in. 

value of x at extreme left end of beam (lower boundary), in. 
value of x at extreme right end of beam (upper boundary), in. 
increment in recurrence solution, * n +± - x n , in. 
rotary inertia, lb-sec 2 
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P matrix (see eg.. (lO)) 

8p matrix (see eq. ( 19 )) 

82 matrix (see eq. ( 20 )) 

£ deflection, in. 

A matrix (see eq. (21)) 

A approximation of second derivative of variable (see eqs. (12) 

to ( 15 ) and. (2l)) 

\|r cross-sectional rotation due to bending, rad 

cd circular frequency, rad/ sec 

increment of cd in trial solutions 


Subscripts: 


A 

"A" solution 

a 

lower boundary 

B 

"B” solution 

l 

upper boundary 

n 

body station 

r 

mode 

Matrix notations: 

0 

column matrix 

1 — 1 
1 1 

square dr rectangular matrix 

-1 

inverse matrix 

L 1 ! 

unit matrix 


Primed symbols denote differentiation with respect to x. 
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DERIVATION OF THE METHOD 


A derivation is presented for the differential equation of an oscillating 
team structure expressed as a set of four first-order differential equations. 

A recurrence formula is established which relates the system variables (mode 
shapes, mode slopes, mode moments, and shears, and their derivatives) at one 
station of the beam to those at an adjoining station. Boundary equations are 
established for the free-free system and a characteristic equation is obtained 
for computing the critical frequencies of the system. Treatment of discontinu- 
ities is considered in detail, and the orthogonality and equilibrium relations 
are established. 


Equations of Motion 

A sketch is presented which shows the equations of motion. 



7 



Consider the differential element of a heam of dx length in the sketch, 
then the summation of the vertical forces on the element will yield 


m 5-L dx + — dx = 0 
dt 2 ^ 


( 1 ) 


Summing moments about p and dropping second-order terms gives 


Vdx + Z^-|dx-^dx = 0 

at 2 Sx 


( 2 ) 


From elementary beam theory, it is known that 


V = 



(3) 


M = El — 

Sx 


(4) 


where t is the cross-section rotation due to bending which differs from 
c^/dx because of shear deformation. 

Also, if interest is confined to the undamped rth natural mode, solutions 
of £, i|f, V, and M are assumed to have the form 


t(x,t) = tr(x)e ia¥t 
V(x,t) = V r (x)e 1C%,t 
M(x,t) = M r (x)e^ a 4'^ 

where is the angular frequency of undamped harmonic motion in the rth 
mode. 


(5) 


It should also be noted that KAG and El are functions of x only; 
therefore 

d(KAG) = d(KAG) 
dx dx 
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c>(El) = d(El) 

ax ^ 

Substituting equations (5) in equations (l) to (4) and performing the indi- 
cated differentiations gives the basic form of the differential equations for 
the natural vibration of a beam, consideration being given to flexure, shear 
deformation, and rotary inertia. Henceforth the subscripts will be dropped for 
purposes of clarity, and the equations will appear as 


-mo)?^ + 

dV 

-r- - ° 
dx 

(6) 

V - 

. ® . 0 

dx 

(7) 

*L + JL 

- Hr = 0 

(8) 

dx KAG 

M d* _ Q 

El dx 

(9) 


Expressed in matrix notation, equations (6) to ( 9 ) appear as 


0 

-mix? 

0 

0 

1 

0 

0 

0“ 


V 

r 



0 

0 

0 

-Zo£ 

0 

1 

-1 

0 




1 

0 

0 

-1 

0 

1 

KAG 

0 

0 


* 

V' 

V 

> = 0 

do) 

0 

0 

-1 

0 

0 

0 

0 

1 

El 


M' 

l M J 




where the primes denote differentiation with respect to x. 

If the 4x8 matrix in equation (10) is denoted as (3, and the 8x1 
column matrix of unknowns as X, then equation (10) may be written as 

ra® ■ ° (ii) 

The matrix equation (eq. (ll)) constitutes the equations of motion of the 
system. 
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Recurrence Solution 


The following relationships are introduced and are based on the assumption 
that the first derivatives of the system variables are linear over the small 
increment x n to x n+ 2.* I*et ( x n4-l “ x n) = Ax- 


d £n+l _ d £n 
dx dx 


+ AxA n (0 




Wl = ST ^ + + Y>n(5) 


( 12 ) 


d t n +l = <^n 
dx dx 


atn 


■\ 


+ AxA n (t) 


*n + l = & + fa + £|t A n W 


( 13 ) 


dV n+ l dVn ^ 

—T— = 7 — + ^n(V) 

dx dx 


V n+ ! = ^ Ax + V n + ^ A n (V) 

J 


( 14 ) 


dMn-i-i dM n , . 

+ AxX n M 

dx dx 




M n +1 = “ ■ ^ + M n + ^ A n (M) 


( 15 ) 


At this point M5), A n (^), et cetera^ are to be regarded as undetermined- 


Equations (12) to ( 15 ) may be written in matrix notation: 
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and 



( 20 ) 


( 21 ) 


Then, from the relationships (17) to (2l), equation ( 16 ) becomes 



The system variables must now satisfy the equations of motion at any specified 
station. That is, when x = x n +q. 



(23) 
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I III III ill I llll III III 


by virtue of the constraint imposed through the following determination of the 
A matrix. Substituting equation (22) in equation ( 23 ) gives 




+ Ax 




= 0 


Solving equation (24) for A n yields 


(24) 





(25) 


Thus, the variables A(£), A(iO, et cetera which were undetermined in equa- 
tions (12) to ( 15 ) are defined by requiring satisfaction of equation ( 23 ). 
Substituting equation ( 25 ) into equation (22) gives the recurrence formula for 
the functions £ , \|r, V, and M: 



It is expedient to perform the matrix manipulations to express equa- 
tion ( 26 ) in the form 


where 



( 27 ) 


(28) 


From equations (id), (ll), and (20) it can be seen that 


HH 


HW* f 


- W* f 


2 


■-1 


For convenience, let 


a = m n+1 o^ 
b = Z n+ iaj2 


d = KAG 


e = El 


n+1 

-1 

n+1 




k - ^ 
151 2 


Then, equation (29) may be expressed as 


HH = 


-ak 

0 

1 

0 


0 

-bk 

-k 

-1 


2 


Ax 


KAG n+1 2 


1 

k 

dk 

0 


0 

-1 

0 

ek 


1 Ax 


EI n +l 2 


The determinant D of this matrix is given by the formula 

D = -aek^ + (l + bek 2 ) (l + adk 2 ) 


( 29 ) 


(30) 


(31) 


(32) 
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By utilizing equations (10), (ll), (19), (20), and (30 ) and by performing 
the matrix operations indicated in equation ( 26 ), tbe elements of the matrix 
p 1 , which appears in equation ( 27 ), may be reduced to the explicit algebraic 

expressions given in equation ( 33 )* Th e elements are functions only of the 
structural properties at x n and x^^, of the interval length Ax, and of the 

unknown-frequency cjd. Equations ( 32 ) and (33) provide the working relationships 
for practical use of the recurrence formula given by equation ( 27 ). 


al^lek^l-bdj-d) 

ak [ek?( 1-bd) -d] 

k 

1 

[ek 2 (l-bd)-djk 

(ek^Cl-bdJ-dj 

ek 2 

ek 

k(l+bek 2 ) 

1+bek 2 

k 2 

k 

{ek 2 ( 1-bd) -dj k 2 

|ek 2 (l-bd)-d]k 

ek 5 

ek 2 

aek 5 

aek 2 

ek 2 [ak?(l-bd)-'0 

ekfal^d-bdj-b] 

ek 2 

ek 

ek(l+adk 2 ) 

e (l+adk 2 ) 

aek 1 * 

aek 5 

k( 1+adk 2 ) 

1+adk 2 

ek 5 

ek 2 

ek 2 (l+adk 2 ) 

ek(l+adk 2 ) 

a( 1+bek 2 ) k 

a( 1+bek 2 ) 

ak 2 

ak 

ak 2 jek 2 (l-bd)-d) 

ak jek 2 (l-bd)-dj 

aek 5 

aek 2 

a(l+bek 2 )k 2 

a ( 1+bek 2 ) k 

ak 5 

ak 2 

k (l+bek 2 ) 

1+bek 2 

aek 11 

aek 5 

ak 2 

ak 

(a^d-bdj-bjk 

(ak^l-bdj-b) 

k 

1 

el^ja^d-bd)-^ 

ek(ak 2 { 1-bd) - 

ak 5 

ak 2 

[I^(l-bd)-b]^ 

(a^d-bdj-bjk 

k 2 

k 

k( 1+adk 2 ) 

1+adk 2 


Boundary Conditions 

Xn the work that follows, the solution will he restricted to considera- 
tions of the free-free vibration behavior. This mode of behavior is appropriate 
to the in-flight characteristics of launch vehicles. Henceforth, the sub- 
scripts a and l indicate that the associated quantity is evaluated at the 
lower boundary x = x a or at the upper boundary x = x^, respectively. The 

boundary conditions at both ends of a free-free beam require that 
Condition (a): 

V a = V* = 0 

Condition (b): 

Ma = M z = 0 
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Substituting condition (a) in equation (8) shows that 


Condition (c): 


Condition (d): 






Substituting condition (a) in equation ( 7 ) gives 


Condition (e): 


Condition (f): 


< = -Z a ^t a 


M { = -Zi^i 


Substituting condition (b) in equation ( 9 ) yields 


Condition (g): 


♦a * = 0 


Solving equation (6) for dV/dx gives 


Condition (h): 


Condition (i): 




Vj » 


These conditions may now be expressed in matrix form. For the lower 
boundary conditions it follows that 
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r c\ 

si 


r ♦ i 

T a 


1 

0 



£a 


0 

1 

< 


0 


0 

0 


>■< 

♦a 

= 

1 

0 


m a^a 

f 

0 


v a 


0 


0 

0 

M a 


-z a o^i a 


-z a a£ 

0 

M a 


0 


0 

0 

V J 


^ J 


_ 

— 



and for the upper boundary conditions 


r > 

s; 


r t, 1 


1 

0 

u 


£ 1 


0 

1 

♦1 


0 


0 

0 


s = < 

* l 

\ = 

1 

0 

Vi 

> 

1 

/ 

0 

m z o£ 

v z 


0 


0 

0 

M i 


-Z 1^1 


-z z a£ 

0 

< > 


0 

k. J 


0 

0 



(34) 


(35) 


Calculation of Frequencies 

Suppose now that a number of stations have been established at definite 
points x n along the span of the beam. The system variables at all stations 

are then uniquely determined if values are assigned to £ a , and 0). 
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Equation (jk) determines the system variables at x = x a . Equation (27) yields 
the system variables at x = x n in terms of the values at x = x a . Once the 
variables at x = x a are known, equation (27) may be used again to compute the 
variables at x = x a +l> an(i ^hls process may be continued until the system 

variables at all stations are known. This process will be termed "a solution 
of the system” for the given values of £ a , t a , and cd. The principle of 

superposition of solutions is employed by using an arbitrary value of o> in 
two separate solutions. The first solution has £ a = ~ 0 and i s called 

the "A” solution j the second has £ a = 0 , ^ a = 1 and is called the "B M solu- 
tion. Then by using conditions (a) and (g) of the section on boundary condi- 
tions, it can be concluded that 



The subscripts A and B designate the separate solutions for the afore- 
mentioned boundary constraints. The coefficients Aq, Ag, Bq, and Bg are 

essentially influence coefficients where Aq and A2 are the values of V 
and V, respectively, at x^ due to a unit value of t at x a . The values 
of V and V, respectively, at x^ due to a unit value of £ at x a are 
termed Bq and B2. 

Nontrivial solutions of equation (36) exist if, and only if, the determi- 
nant of the square matrix is zero. The natural frequencies are therefore deter- 
mined by the equation 


Al 

Bl 

a 2 

b 2 


( 37 ) 


hereinafter referred to as the frequency equation f ( cd) . The coefficients of 
equation (37) are frequency dependent as is evident from inspection of the 
P n +i matrix defined by equations (30), ( 32 ), and ( 33 )- This frequency depend- 
ency permits the satisfaction of equation (37) and yields the natural frequen- 
cies of the system. Also, from superposition of the "A" and "B tT solutions, the 
following additional relationships can be stated: 
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I 



( 38 ) 



^Z A 

^ZB 

+ZA 

^ZB 



Cl D x 


Dr- 



where C]_, C 2 > I>|_, and Dg are influence coefficients for and \|fj , 

From the second row of equation (36) 


A 2^a + B 2^a = 0 


(39) 


Hence, 


*, - - 


Ms 

A 2 


(40) 


From the first row of equation ( 38 ) 


U = C^a + D^ a 

Solving equation (4l) for f^ a and substituting equation (40) for >|r a 


(41) 


gives 


a 2^ 2 

- B2C ]j_ 


and normalizing in terms of = 1 yields 


t = a 2 
ba a 2 d 1 - b 2 c 1 

Substituting equation (42) in equation (40) yields 


(42) 


♦a 


-B 2 

AgDi “ B2C1 


(43) 


It is reminded that the appropriate values of Ap, A 2 , Bp, B 2 , Cp, 

C 2) Dp, and D 2 are only those associated with "A" and "B" solutions 
satisfying equation (57)* The numerical procedure in obtaining valid solutions 
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is to perform trial solutions for assumed values of a) until equation (37) is 
satisfied. The lowest frequency satisfying the equation yields the fundamental 
natural mode of vibration. Succeeding higher frequencies, which satisfy equa- 
tion (37) , yield the progression of overtones or higher modes. 


Mode Shapes and Related Data 

The coefficients to Dg associated with the proper values of cd 

which satisfy equation (37) > when substituted into equations (42) and (43), 
yield appropriate initial values of £ a and i|f a for the mode shapes. The 
recurrence formula of equation ( 27 ) relates the modal characteristics at a given 
station to an adjacent station. Thus, having the initial values leads to a 
complete solution for all stations. Such a solution yields the mode shapes, 
cross-section rotations, mode moments, mode shears, and the first derivatives 
with respect to x of each function. 

It should be noted that in letting ^ equal unity in obtaining equa- 
tions (42) and (43), the system variables derived from use of the two relation- 
ships are henceforth compatible with the mode shape normalized at x = xj 
(i.e., £2 = 1 ). For purposes of preventing conflicts in units, it should be 

understood, that the form of the system variables, when = 1, is in essence 



(MO 


That is, £, V, and M from a final solution are actually the ratios 

^>nAz> ^nAl > ^nAz* an( i ^nAz> respectively. These ratios are the system 
variables per -unit of mode deflection at x = xj . 


Treatment of Discontinuities 

In most practical problems associated with nonuniform beams, El, KAG, 
Z, and m exhibit many discontinuities over the range x a x ^ X 2 • Conse- 
quently, provisions must be made for a convenient solution across discontinu- 
ities. When a discontinuity is crossed, the system variables have to be 
reevaluated. 

From physical consideration M, and V are continuous, and it 

will be shown that \|r’ , £ r , M* , and V f are discontinuous. 

Let the subscripts (-) and ( + ) denote values of system variables at a 
discontinuity which is approached through values of x lower than or higher 
than, respectively, the value of x at the discontinuity. 
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Then from equation (9) where M( _ j = M(+) 


dt 


(+) . M (-) 


dx EI( + ) 

From equation (8) where = V(+) and ^(-) = 


*(+) 


(^ 5 ) 


and from equation (7) 


45 (+> 


dx 


= t ( _) - 


(-) 


KAG 


(+) 


dM 

dx 


l±i = T ( _) - Z ( + ) ® 2 t(.) 


( 46 ) 


( 47 ) 


From equation (6) with ^ 


dV 


( + ) _ 


dx 


= »( +) ^S ( .) 


(M) 


The relationship across a discontinuity may then he shown in matrix form 



L° 


o 


m / , \o£ 


w 

0 

0 

0 


0 


0 

0 

0 

0 

0 


1 

0 

0 

' J (+) c 

0 


-Li . \o£ 


0 


0 


0 

0 

0 

0 

0 


KAG(+) 
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0 

0 

0 

1 

1 


0 

0 

0 

0 

0 

0 

0 


0 

1 


EI (+) 

0 

0 

0 

0 



(^9) 
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This relationship is exactly satisfied when Ax = 0 is entered into P n+ i 

of equation ( 33 )* This contributes to the efficiency of the method for use in 
conjunction with high-speed digital computers because the same recurrence for- 
mulas are used to relate system variables across a discontinuity as are used to 
relate variables across a finite interval. 


Orthogonality Condition 

One of the major uses of natural mode data is in modal form series solu- 
tions. Their practical applications in series solutions is a result of the 
many simplifications that result from their orthogonality relationships. These 
relationships for modal solutions, restricted to elementary beam theory, are 
well known and widely referenced in reports and standard texts on elementary 
vibrations. The addition of rotary inertia and shear deformation alters the 
conventional relationship and a clear description of the orthogonality rela- 
tionship for nonuniform free-free beams with secondary effects is considered 
essential to the completeness of this paper. In reference 6 , Leonard derives 
the orthogonality relationships for the uniform free-free beam with shear 
deformation and rotary inertia. In the following discussion, the conditions 
for the nonuniform free-free beam will be set forth. 

Substituting equation ( 8 ) into equation ( 6 ) and adding subscripts to £ 
and \|r to designate the specific mode r yields 




dx 


KA.G 


f <%x 

idx 


= 0 


( 50 ) 


Also, substituting equations ( 8 ) and (9) in equation (7) gives 




r 1 2 ,u cl L-r d ^r 


. Zo^ r - ^EI —J » 0 

Multiplying equation (50) by £ n and (51) by ^ n , adding both results 
together, and integrating over the length of the beam, yields 


(51) 




r x i / 

\ r x z 

d 


'Hr \| 

/ Kn^r + dx = - / 

& 

kag| 

.ST - 

x a \ 





dx 


- 1 T s( ei s 1 )* 1 

a 

r x z ( d£ r 

- J KA&Un ‘ trtn 

x a V > 


dx 


(52) 
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Since the r and n subscripts have been selected arbitrarily, equa- 
tion (52) is equally valid when expressed by reversing r and n. Making 
this reversal yields 



Expanding equations ( 52 ) and (53) by integrating, by parts, the first and 
second terms on the right-hand side of the equations and then subtracting the 
extended version of (53) from the extended form of ( 52 ) will yield 



For the free-free boundary conditions 


kag(— ^ - 

ydx 


El 


gjn 

dx 


x l 


^a 



= V n (xa) - V n (x 2 ) = 0 


Equally valid 
for n or r 


= M n( x i) - M n( x a) = 0 


J 


(55) 


With these relationships, equation (54) leads to the final conclusion that 





") f = 0 

x a 


(56) 
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and when ^ it is necessary that 



< n ?r + z Vr) to = 0 


<% / <%• ( 57 ) 


When it is usual to define the preceding integral as the effec- 

tive mass of the mode. 


% = J Z ^ r 2 + Z\|r r 2 jdx = (br (58) 

To avoid confusion and to obtain the effective mass in terms of mass units 
(lb-sec 2 /in . ), it is necessary only to interpret the symbols and i|f r as 

applicable to the dimensionless mode. 

Now, instead of subtracting as was done to obtain equation ( 54 ), the equa- 
tions are added; then with the additional results of equation (57) when n / r, 
it is found that 



«!> + raol^S 

dx dx Vdx 




n ± r (59) 


and when n = r, it is usual to define the preceding integral as the effective 
spring constant of the mode: 



d ^ r 2 

dx 


+ KAG 




n = r (60) 


where the same statement regarding units for equation (58) applies. 

Equations ( 57 ) and ( 59 ) ar ® the orthogonality properties of the free -free 
beam of variable mass and stiffness with shear deformation and rotary inertia. 

Equations (58) and (60) are useful relationships defining the effective 
mass and effective spring of the rth mode referenced at the normalizing point 
of the mode shape £ r = 1. 
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II 


Equilibrium Relationships 

In addition to the useful orthogonality relationships, two other valuable 
integrals associated with the equilibrium of the mode are frequently encoun- 
tered. For the free-free system to possess dynamic equilibrium while under- 
going natural vibrations, the following relationships prevail: 



m£ r dx = 0 


(61) 




(62) 


SOME ASPECTS OF APPLICATIONS 


Basic relationships for the recurrence solution of the beam vibration 
problem are given in equations (27), (30), (32), and ( 33 )- The natural fre- 
quencies are obtainable from equation ( 37 )* Orthogonality and equilibrium 
relationships are given in equations (57), ( 59 ) ^ (6l), and (62). 


Descriptiveness of Output Data 

The type of solution set forth by the foregoing equations is inherently 
adaptable to producing highly descriptive modal data. It has been shown by 
equation (^ 9 ) that f 1 , £*, M* , and V f are discontinuous functions by vir- 

tue of discontinuities in El, KAG, Z, and m. The quantities obtained by 
the solution are of equal descriptiveness to that of the input data; that is, 
for every discontinuity in input data (El, KA.G, Z, and m) there will result 
an associated discontinuity in the output functions. The quality of these data 
cannot be equaled by the discrete mass methods of analysis such as that of ref- 
erence 3- The method is also inherently suitable for automatic digital- 
plotting routines since the input can be supplied or generated by interpolation 
to a definition as fine as desired to give an acceptable plot of output. 


Numerical Solutions 

The successful use of the foregoing formulation is largely dependent on 
the details of the computer program. Several important features to be observed 
are set forth in the following text. 

Machine time .- Improper programing can result in lengthy computing times 
yielding unreasonably costly data. Conversely, proper programing for minimum 
computer time can be most rewarding with the subject formulation. For example, 


1 
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programing the inverted and extended form of equation (28) given by equa- 
tions (32) and (33) reduced the machine time to one-fourth that required by a 
program performing internally the operations of equation (28). Unduly stringent 
accuracy requirements in determining the natural frequencies by equation (37) 
can also result in exorbitant machine time. Extensive use of the method of 
this paper has indicated that three modes with frequency convergence to 0.1 per- 
cent error can be achieved with around 300 stations of input data to an IBM 7094 
electronic data processing system in about 5 minutes. It should be noted that 
this accuracy on frequency is valid for the mathematical model but it is not 
indicative of the accuracy of predicting the frequencies of the actual struc- 
ture. The accuracy is obviously dependent on the appropriateness of the math 
model. 

Attempting to satisfy upper boundary conditions to unneccessary accuracies 
can also be costly. For example, \|r* , Vj, and are zero for the free-free 

boundary case. These absolute conditions cannot be achieved by digital com- 
puters. The analyst, however, should accept finite boundary values that are a 
fraction of 1 percent of the peak absolute values of their respective functions 
over the total span x a < x < x^ . Boundary values which are zero correct to 

the fifth significant figure of the maximum value of a function can readily be 
obtained. 

Near zero stiffness .- In applying the recurrence technique of this report 
the unwary analyst can be confused and misled by using inputs of KAG and El 
that are near zero at or near the upper boundary position. This condition is 
encountered particularly with pointed nose cones or other similar structures 
that essentially taper to zero. The problem manifests itself in a radical 
variation of the mode shape in the area approaching the upper boundary of the 
span. This variation is a combined result of near zero stiffness of KAG 
and/or El and the failure to achieve numerically the absolute theoretical 
zero boundary condition. This condition is readily corrected by avoiding, near 
the free ends, KAG and El input values of less than 0.0001 of their respec- 
tive average values over the total span. 

Supe rposition for final modes .- The sequential nature of a solution by 
the method proposed in this paper is as follows: First, obtain for a variety 

of frequencies the influence coefficients Aq, Ag, Bq, and B2 by the ft A !t 
and "B" solutions discussed in "Calculation of Frequencies." From this trial 
technique, a frequency is discerned to a desirable accuracy that will satisfy 
the frequency equation (eq. ( 37 ))- Second, after obtaining a natural frequency, 
the influence coefficients associated with the critical frequency are substi- 
tuted into equations ( 42 ) and ( 43 ) to obtain the boundary values £ a and 

associated with the normalized natural mode. Finally, with the knowledge of 
an(i ^a> initial conditions and coefficients can be fully ascer- 

tained, and the recurrence solution of equation (27) can be extended over the 
total span. This solution yields data on all of the system variables for all 
desired stations. An important feature associated with this last operation 
should be noted. Entering £ a and i|r a simultaneously into equation (27) in 

the final solution for the modal characteristics is a different digital opera- 
tion than the original procedure of superposition that led to knowledge of £ a 


27 



and ty a . Consequently, a simultaneous solution will result in different dig- 
ital round-off errors that will reduce the accuracy of the final boundary 
values below that previously established by superposition. Since a large num- 
ber of significant figures are essential to the mechanism of the numerical 
solution, a duplicate procedure for the initial operation is recommended in 
obtaining the final modal data. The available values of £ a and should 

be used in independent "A” and T, B Tt solutions and should be combined at corre- 
sponding x stations to obtain the desired final modal data. 

Integration intervals . - The fundamental differential equations of motion 
given by equation (10) are all of first order. Integration is performed on all 
variables with the assumption of linear variations of the first derivatives as 
defined by equations (12) to ( 15 ). For this assumption to be acceptable, the 
Ax interval must at all stations be suitably small. Studies on interval sizes 
have shown that adequate results are normally obtainable by describing the 
physical characteristics of a vehicle at all points of discontinuity with the 
added constraint that Ax = L/lOO. The latter constraint is easily programed 
into the computer and obviates a large amount of repetitive input data over 
long constant sections. In addition, where continuous but nevertheless radical 
variations occur in the input functions, additional stations should be included 
to insure that the variations are adequately defined. 

Frequency search technique .- The economy of the machine solution is inti- 
mately associated with the iterative technique employed to obtain the unique 
frequencies that will satisfy equation (37)* Equation (37) is of the form 

f (co) = 0 

and must be solved by trial and error. 

Successful results have been achieved by performing repeated trial solu- 
tions for f (co) by systematically increasing the frequency 0 ) by AjD until a 
sign change is indicated, i.e.. 


f(aP) 

where aJ 1 is the nth trial value. At this point, a second-order (parabolic) 

curve fit is applied to points at f^opj, and and the 

resulting analytical expression can be solved for co where f(co) = 0. This 
operation can be continued with each improved frequency until the difference 
between successive frequencies is less than some preassigned tolerance. 

Convergence can be hurried by good estimates of An) for the specific 
problem and refraining from a too stringent tolerance in frequency. Increasing 
An) and the frequency tolerance for successive overtones is also advantageous 
in maintaining machine time within practical limits. 
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Secondary Influences 


A principal objective of the mathematical development of this paper was to 
provide a means for including, as well as investigating, the effects of the sec- 
ondary contributions of shear deformation and rotary inertia in modal analyses. 

Comparison with classical cases .- Studies have been made by comparing the 
classical or exact modal data for uniform free-free beams without secondary 
effects with data obtained by the subject recurrence technique with secondary 
effects. The two methods, secondary effects being excluded, agreed to an accu- 
racy of greater than 0.01 percent on the first three mode frequencies, mode 
shapes, and nodal points. It was felt that this agreement qualified the numer- 
ical technique for use in further comparative studies. 

The salient results of some studies to determine the significance of shear 
and rotary inertia on the classic uniform beam are furnished in figure 1. Data 
were derived for a solid cylin- 
drical beam and for a thin-wall 
cylindrical beam of a thickness- 
to-diameter ratio of 0 . 015 . 

Since the effectiveness of shear 
deformation is dependent upon 
the ratio of the shear modulus 
G to the bending modulus E, 
the studies were made for both 
fiber-glass and steel beams. It 
can be seen from the curves for 
the steel beam, that for either 
of the cross sections the com- 
bined secondary influences 
caused frequency reductions in 
the first three modes of not 
greater than 5 percent when 
ratios of length to diameter 
(L/D) are greater than 20. This 
frequency reduction compares 
with those reductions of up to 
Ik percent for values of L/D 
greater than 20 for fiber -glass 
beams. The value of G/E for 
the fiber-glass beam was 0.20 
whereas that for the steel beam 
was 0.^2. It is also evident 
that the frequency reductions 
for both the steel and fiber- 
glass beams are more pronounced 
for the thin-wall cases than for 
the solid beams. As the aspect 

ratio (L/D) decreases, the reduction in frequency becomes quite pronounced and 
reaches more than 50 percent in the third mode for the thin-wall cases for 
aspect ratios of k- and 7 for the steel and fiber-glass beams, respectively. 


Length-to-diameter ratio, L/D 



(a) Steel beam. 

Length-to-diameter ratio, L/D 



(b) Fiber-glass beam. 

Figure 1.- Effects of rotary inertia and 
shear deformation on the frequencies 
of free-free vibrations of uniform 
cylindrical beams. 
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In reference 7, Kruszewski presented closed-form solutions for the uni- 
form beam with shear deformation and rotary inertia. The report gives general 
data for readily determining the significances of the secondary influences for 
general beam cross section and elastic properties. The curves of figure 1 
obtained by the numerical procedure of this paper were checked against the data 
from the closed-form solution of reference 7 and were found to be in nearly 
perfect agreement. This comparison further qualifies the numerical technique 
as a satisfactory and accurate procedure for incorporating the secondary influ- 
ences into the beam solution. 

Effects on launch vehicle .- The six data points given in the fiber-glass 
curves of figure 1 show the frequency reductions due to secondary influences on 
a four-stage space vehicle with L/D = 25 and the fourth stage of a space 
vehicle with L/D = 5* Of the first-stage length 5*5 percent was fiber glass 
and 2 6 percent of the fourth-stage length was fiber glass. The first-stage data 
reveal reductions greater than those predicted for thin-wall steel beams but 
are less, as would be expected, than those for thin-wall fiber-glass beams. 

The fourth-stage frequencies for L/D = 5 were significantly reduced by 
the secondary effects and closely correspond to the reductions predicted for the 
thin-wall fiber-glass beam. It appears that the fiber glass appreciably con- 
tributed to the frequency reduction even though the fiber glass extended over 
only one -fourth of the stage length. The fiber-glass section extended over a 
span running from 0.25 to 0.55 of the length which was a region subjected to a 
significant proportion of the shearing action. The specific vehicle for which 
these data were derived is documented further in the section entitled 
"NUMERICAL EXAMPLE." 

Measure of significance .- The aforementioned considerations suggest the 
usefulness of the curves of figure 1 in estimating the probable significance of 
shear deformation and rotary inertia in free-free natural frequencies of launch 
vehicles. Both the vehicle data and the curves clearly indicate the necessity 
of including secondary influences if accurate analyses are desired on low-aspect 
ratio s tinctures. 


NUMERICAL EXAMPLE 


An example is presented of an application of 
the actual research vehicle illustrated in figure 


First-stage separation Second-stage separation 


Third-stage separation 



Figure 2.- Typical multistage research 
launch vehicle. 


the recurrence solution to 
2. Solutions have been pro- 
vided for the first- and 
fourth-stage configurations 
of flight. These two cases 
are submitted to show an 
actual numerical example 
derived from the subject 
method and to illustrate 
the varying importance of 
secondary effects with the 
vehicle’s length-to-diameter 
ratio. The basic physical 
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characteristics of the vehicle required as input to the program have been 
recorded and submitted in table I* The tabulation provides input for both 
cases by taking the physical characteristics at x = 619*85 as the first 
applicable quantities for the fourth-stage analysis. 


Frequency Data 

In table II , the frequencies and percent reduction in frequencies are 
given for the first three modes of vibration for the two stages of the vehicle 
of figure 2 . The data are displayed for analyses without the secondary effects, 
with both effects, with rotary inertia only, and with shear deformation only. 

Inspection of the data for the first stage will show only small effects of 
shear deformation and rotary inertia. The influences grow progressively with 
increasing modes as would be normally anticipated. The maximum contribution of 
rotary inertia and shear deformation is seen to be a reduction in frequency of 
slightly less than 5 percent on the third mode. The first-stage configuration 
has a length-to-diameter ratio of approximately 25 - In the fourth- stage con- 
figuration with an aspect ratio of 5 - 0 , the importance of rotary inertia and 
shear deformation are quite evident. The combined secondary influences produce 
a reduction of 58 percent in values computed without consideration of shear 
deformation and rotary inertia. The significance also increases with increasing 
modes. These results strongly suggest that for reliable frequency calculations 
on the upper stages of typical launch vehicles, consideration of the secondary 
contributions to flexure must be given. 

It is interesting to note, however, that good results would have been 
obtained in all cases shown by incorporating in the solution the effects of 
shear deformation only. The frequency error most affected by the secondary 
influence (third mode of the fourth stage) would be slightly less than 7 Per- 
cent if shear deformation only had been considered. 


Elastic Curve Characteristics 

A graphical comparison of the modal functions for the first three natural 
modes for the first stage of the example vehicle would reveal only trivial 
departure from concurrency. However, it will be shown later that this is not 
the case for the fourth-stage configuration where the secondary effects have 
been shown to be appreciable. 

Mode shapes .- In figure 3 * the fourth-stage mode shapes £ are compared 
for solutions with and without secondary effects. The comparison of mode 
shapes was accomplished by the application of the method of least squares. In 
spite of the severe reductions noted in frequencies, the first three mode 
shapes retain similar graphical characteristics. 

Mode shape inflections near tips .- An interesting elastic curve character- 
istic should be noted that results from the inclusion of shear deformation. 

For beams with secondary effects, an inflection in the mode shape is frequently 
observed near one or both of the free ends and outboard of the outermost nodal 
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Figure 3-- Effects of rotary inertia 
and shear deformation on mode 
shapes of a low-aspect-ratio 
upper stage of a launch vehicle. 


point. Inspection of either of the 
three modes of figure 3 will indicate 
the inflection or reversal in curva- 
ture on the right-hand end in the 
vicinity of x = 90 inches. The 
inflection point occurs when the sec- 
ond derivative of the mode shape 
becomes zero. An expression is 
readily obtained by differentiating 
equation (8) that indicates the 
parameters involved in the phenomenon, 
that is 


since 


d 2 £ d/ V \ dt _ 
^2 + dx\KAG J " dx “ 



(63) 


at the inflection point, and 



< 0 


to the right of the inflection point. Then it follows from equation (63) that 
for the dip in the elastic curve on the end to occur 


A. JL > 2 ( 6k ) 

dx KAG dx 

Equation (64) simply states that when curvature contribution due to shear 
deformation is greater than the curvature contribution to the elastic curve due 
to bending, the inflection near the tip will be evident. In elementary beam 
solutions, shear deformation is not considered) this case corresponds to the 
limiting case where KAG -» <». It is then obvious from equation (6k) that for 
cases where shear deformation is ignored, the shear deformation contribution to 
curvature is zero and the conditions of equation (64) can never be realized. 

Mode slopes . - A pronounced departure between data computed with and with- 
out secondary effects is seen in the comparison of mode slopes in figure 4. 

The solid curves are the slopes computed without shear deformation and rotary 
inertia. The dashed curves were computed with secondary effects. The slopes 
are compatible with the amplitudes of the mode shapes given in figure 3 as 
established by the least-squares method. The most conspicuous differences are 
the discontinuities that are evident in the data incorporating rotary inertia 
and shear deformation. The discontinuities result from shear deformation and 
are not evident in results without shear deformation, whether Including rotary 
inertia or not. Accurate mode slopes frequently are required in system 
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stability studies of autopilot 
structural feedback. It is evident 
from figure k that, for accurate 
structural feedback analyses on 
short L/D configurations, the sec- 
ondary effects, particularly shear 
deformation, should be Included in 
developing the required mode slopes. 

Mode slopes compared with 
cross-section rotation .- In elemen- 
tary beam theory, which ignores both 
rotary inertia and shear deforma- 
tion, the cross-section rotation is 
quantitatively the same as the 
slope. The differences between the 
slope of the elastic curve £ 1 and 
the cross-section rotation are a 
further measure of the significance 
of secondary effects. A comparison 
of these modal functions is shown 
graphically in figure 5 for the 
fourth stage of the example vehicle. 
The observed differences between 
£ 1 and ^ are completely a 
result of shear deformation and 
are not evident in solutions 
ignoring shear. 


Load Characteristics 



1 ! ! 1 L 

20 bo 60 80 100 


x, Inches 


Figure 4.- Effects of rotary inertia 
and shear deformation on mode 
slopes of a low-aspect-ratio 
upper stage of a launch vehicle. 



of Modes 


The moment and shear dis- 
tributions for the fourth stage 
of the example vehicle of fig- 
ure 2 are submitted in figures 6 
and 7« Comparative curves are 
given for data calculated with 
and without the secondary influ- 
ences of shear deformation and 
rotary inertia. The moment and 
shear values are compatible with 
the relative mode amplitudes shown 
in figure 3 and as established 
by the least- squares method of comparison. Both the moment and shear data show 
significant differences with increasing modes between solutions with and with- 
out secondary Influences, the magnitudes of the differences increasing with 
increasing modes. An appreciable separation in the right nodal point of the 
third mode moment curves is observed. The startling differences that are 
observed between amplitudes of the comparative data for the second and third 



Figure 5 .- Comparison of cross-section 
rotations and mode slopes of a low- 
aspect -ratio upper stage of a launch 
vehicle. (All curves include the 
effects of shear and rotary 
inertia. ) 
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modes are principally attributed to the lower frequencies and consequently 
lower inertia loading on the beam for the solutions with secondary effects. 



Figure 6.- Comparison of the mode 
moments from solutions with 
and without secondary effects 
for a low-aspect-ratio upper 
stage of a launch vehicle. 



x, inches 

Figure 7*- Comparison of mode shears 
from solutions with and without 
secondary effects for a low- 
aspect -ratio upper stage of a 
launch vehicle. 


CONCLUDING REMARKS 


A recurrence solution is presented that is especially suitable for 
obtaining highly descriptive modal data on severely discontinuous nonuniform 
beamlike structures; the solution includes the secondary influences of rotary 
inertia and shear deformation. Data are provided for helping the analyst in 
estimating the probable significance of shear deformation and rotary inertia 
in free-free natural frequencies of launch vehicles. 

Numerical examples of the solution are included for a first-stage and a 
fourth-stage configuration of a research vehicle. The secondary influences of 
rotary inertia and shear deformation had trivial effects on the first-stage 
data, but were found to have very significant effects on the data for the 
fourth stage. 

In most practical applications to launch vehicles, it is concluded that 
shear deformation is considerably more important than rotary inertia. Accept- 
able results could generally be obtained by including shear deformation only. 
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The authors’ experiences indicate that, in general, for most team applica- 
tions the secondary influences have the least obvious effect on the mode shapes. 
An appreciable influence is noted in reducing frequencies and large effects are 
experienced on the mode shears and moments. The appearance of discontinuities 
in the mode slopes is a further influence resulting from the consideration of 
shear deformation. 

The addition of shear deformation to elementary beam vibration theory can 
result in an inflection in the mode shape outboard of the outermost nodal 
points. 

Successful and economical solutions by the subject method are intimately 
dependent upon details of the computing program. Some salient points are dis- 
cussed on appropriateness of input data, superposition problems, numerical 
integration intervals, and frequency iteration routine. 


Langley Research Center, 

National Aeronautics and Space 
Langley Station, Hampton, 


Admini strati on , 

Va., April 9 , 1965. 
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TABLE I.- PHYSICAL CHARACTERISTICS OF A MULTISTAGE LAUNCH VEHICLE 


in. 

Z, 0 
lb-sec 2 

m. 

lb-sec 2 /in . 2 

El, 

lb -in . 2 

KAG, 

lb 

-54.35 

1.70 

0.01743 

0.100 x io 9 

10.00 X 10 6 

-19.55 

1.70 

.01743 

.100 

10.00 

-19.55 

14.70 

.06030 

14.770 

33.5^ 

0 

14.70 

.06030 

14.770 

33.5^ 

17.85 

14.70 

.06030 

6.743 

33.5^ 

17.85 

14.70 

.06030 

6.7^3 

33.5^ 

19.55 

14.70 

.06030 

6.7^3 

33.5^ 

19.55 

15.59 

.13094 

6.743 

59.3^ 

20.85 

15.59 

.13094 

6.7^3 

59.3^ 

20.85 

15.59 

.13094 

3^.59 

59.3^ 

71.65 

15.59 

.13094 

34.59 

59.3^ 

71.65 

6 . 6 l 

.11350 

3^.59 

59.3^ 

204.85 

6 . 6 l 

.11350 

3^.59 

59.3^ 

204.85 

6 . 6 l 

.11350 

8.352 

59.3^ 

206.35 

6.61 

.11350 

8.352 

59.3^ 

206.35 

1.50 

.00895 

8.352 

209.00 

207.85 

1.50 

.00895 

8.352 

209.00 

207.85 

1.50 

.00895 

94.92 

209.00 

209.95 

1.50 

.00895 

94.92 

209.00 

209.95 

1.50 

.00895 

6.914 

209.00 

211.45 

1.50 

.00895 

6.914 

209.00 

211.45 

.55 

.00895 

6.914 

26.46 

212.95 

.5^ 

.00895 

6.914 

26.00 

212.95 

.5^ 

.00895 

13.900 

26.00 

217.95 

.51 

.00895 

10.400 

21.00 

222.95 

.20 

.00895 

8.000 

17.00 

227.95 

.12 

.00895 

6.000 

i 4 . 8 o 

232.95 

.05 

.00895 

3.200 

15.60 

236.35 

.047 

.00895 

2.100 

15.10 

236.35 

.047 

.00895 

1.403 

15.10 

237.85 

.045 

.00895 

1.403 

12.92 

237.85 

1.40 

.01474 

1.403 

31.85 

239.55 

1.595 

.01474 

1.403 

36.89 

259.55 

1.595 

.01474 

4.000 

36.89 

245.25 

1.556 

.01474 

3.370 

56.71 

245.25 

1.556 

.01474 

1.042 

56.71 

246.75 

1.321 

.01474 

1.042 

61.75 

246.75 

1.321 

.01474 

1.042 

61.75 

248.25 

1.316 

.01474 

1.042 

61.75 

248.25 

1.316 

.01474 

3.361 

61.75 

257.85 

1.200 

.01474 

3.361 

61.75 

257.85 

1.200 

.03673 

3.361 

61.75 

258.75 

1.200 

.03673 

3.361 

61.75 

258.75 

1.200 

.03673 

4.910 

33.69 

262.85 

1.200 

.03673 

4.910 

33.69 

262.85 

.585 

. 026668 

4.910 

33.69 

424.85 

.385 

.026668 

4.910 

33.69 

424.85 

.385 

.026668 

1.665 

33.69 

426.35 

.385 

.026668 

1.665 

33.69 

426.35 

1.090 

.011718 

1.665 

24.88 

427.85 

1.090 

.011718 

1.665 

24.88 
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TABLE I.- PHYSICAL CHARACTERISTICS OF A MULTISTAGE LAUNCH VEHICLE - Concluded. 



z. 

m, 

lb-sec 2 / in . 2 

El, 

KAG, 

in. 

lb-sec 2 

lb -In . 2 

lb 

427.85 

1.090 

0.011718 

3.025 x 109 

24.88 x 106 

444.85 

1.090 

.011718 

3.025 

24.88 

444.85 

1.090 

.011718 

1.071 

24.88 

446.55 

1.090 

.011718 

1.071 

24.88 

446.35 

.385 

.03360 

1.071 

33.69 

447.85 

.385 

.03360 

1.071 

33.69 

447.85 

.385 

.03360 

4.910 

33.69 

451.55 

.385 

.03360 

4.910 

33.69 

451.35 

.385 

.026668 

4.910 

33.69 

604.55 

.385 

.026668 

4.910 

33.69 

604.55 

.385 

.026668 

.5224 

33.69 

606.05 

.385 

.026668 

.5224 

33.69 

606.05 

.084 

.00585 

.5224 

24.08 

! 607.55 

.091 

.00585 

.5224 

24.61 

607.55 

.091 

.00585 

3.630 

24.61 

610.00 

.103 

.00585 

4.080 

25.48 

6 l 4.00 

.122 

.00585 

4.810 

26.89 

618.35 

.143 

.00585 

5.600 

28.43 

618.35 

.143 

.00585 

.8931 

28.43 

619.85 

.150 

.00585 

.8931 

28.96 

619.85 

.070 

.00266 

.8931 

6.639 

621.35 

.070 

.00266 

.8931 

6.639 

621.35 

.070 

.00266 

1.409 

6.639 

640.95 

.070 

.00266 

1.409 

6.639 

640.95 

.070 

.00266 

.2714 

6.639 

642.45 

.070 

.00266 

.2714 

6.639 

642.45 

.900 

.04457 

.2714 

1.169 

643-95 

.900 

.04457 

.2714 

1.169 

643.95 

.900 

.04457 

.3760 

1.169 

647.45 

.900 

.04457 

.3760 

1.169 

647.45 

.900 

.04457 

.3760 

1.169 

667.96 

.900 

.04457 

.3760 

1.169 

667.96 

.900 

.04457 

.3760 

1.169 

671.462 

.900 

.04457 

.3760 

1.169 

671.462 

.900 

.04457 

.20014 

1.169 

672.962 

.900 

.04457 

.20014 

1.169 

672.962 

.140 

.00344 

.20014 

1.000 

674.462 

.140 

.00344 

.20014 

1.000 

674.462 

.140 

.00344 

.250 

1.000 

698 . 45 

.140 

.00344 

.250 

1.000 

698.45 

.i4o 

.01409 

.250 

1.000 

699-76 

.140 

.01409 

.250 

1.000 

699.76 

.140 

.01409 

.1364 

1.000 

701.26 

.l4o 

.01409 

.1364 

1.000 

701.26 

.380 

.01409 

.1364 

.500 

702.76 

.380 

.01409 

.1364 

.500 

702.76 

.380 

.01409 

.250 

.500 

712.65 

.380 

.01409 

.250 

.300 

712.65 

.380 

.000259 

.250 

.500 

716.95 

.380 

.000259 

.250 

.500 
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TABLE II.- REDUCTION OF FREE-FREE NATURAL FREQUENCIES OF A MULTISTAGE RESEARCH VEHICLE 


DUE TO ROTARY INERTIA AND SHEAR DEFORMATION 


Frequency, radians/sec 


Reduction in frequency, percent 


Without rotary With rotary With rotary with shea r 


inertia and 


inertia and 


I shear deformation I shear deformation 


inertia only deformation only 


inertia and 
shear deformation 



Four-stage launch vehicle, aspect ratio » 25 




Fourth stage of multistage launch vehicle, aspect ratio « 5.0 



400.968 

348.945 

23.66 

8.57 

20.43 

1104.970 

803.663 

44.38 

19.19 

41.22 

1655.463 

1049.897 

58.09 

27.78 

54.20 
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"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof .” 

— National Aeronautics and Space Act of 1958 
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